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Abstract. Mixing a passive scalar field by stirring can be measured in a 
variety of ways including tracer particle dispersion, via the flux-gradient rela- 
tionship, or by suppression of scalar concentration variations in the presence of 
inhomogeneous sources and sinks. The mixing efficiency or efficacy of a partic- 
ular flow is often expressed in terms of enhanced diffusivity and quantified as 
an effective diffusion coefficient. In this work we compare and contrast several 
notions of effective diffusivity. We thoroughly examine the fundamental case of 
a steady sinusoidal shear flow mixing a scalar sustained by a steady sinusoidal 
source-sink distribution to explore apparent quantitative inconsistencies among 
the measures. Ultimately the conflicts are attributed to the noncommutative 
asymptotic limits of large Peclet number and large length-scale separation. We 
then propose another approach, a generalization of Batchelor's 1949 theory of 
diffusion in homogeneous turbulence, that helps unify the particle dispersion 
and concentration variance suppression measures. 

1. Introduction. Flow-enhanced mixing is an important phenomenon in natural 
systems varying in size from as small as human cells to as large as the atmosphere 
and the ocean and beyond [5, 6, 15]. The enhancement of molecular mixing by 
stirring can be observed even for simple laminar flows, and a quantitative under- 
standing of fundamental mechanisms and properties of mixing processes is key to 
accurate modeling of these systems. 

Passive scalars are mathematically idealized entities that serve as proxies to 
formulate and investigate this problem. Given its initial location, the trajectory 

2000 Mathematics Subject Classification. 76F25, 76M50, 74Q20, 76R05. 
Key words and phrases. Transport &; mixing; Stirring & mixing; Turbulent diffusion. 
This work was supported in part by NSF awards DMS-0553487, PHY-555324, and PHY- 
0855335. 



2 



ZHI LIN, KATARINA BODOVA AND CHARLES R. DOERING 



of a passive tracer in R d advected by a prescribed flow field u(x, t) is described by 
the stochastic differential equation 

dX(t) = u(X(t),t)dt + dW(t), (1) 

where k is the molecular diffusivity and W is canonical d-dimensional Brownian 
motion. Equivalently, the Fokker-Planck equation that governs the evolution of the 
scalar concentration field T(x, t) stirred by an incompressible (V • u = 0) flow is the 
advection-diffusion equation 

d t T + u • VT = kAT (2) 

supplemented with initial concentration T(x, 0) and appropriate boundary condi- 
tions. In many applications, the scalar field is constantly replenished and depleted 
by sources and sinks. Consequently, the homogeneous partial differential equa- 
tion (2) would be appended with an inhomogeneous term corresponding to the 
source-sink distribution whose relationship with the stirring further adds to the 
mathematical complexity of the problem. 

In this work we first review several different measures established in the literature 
to characterize the effects of stirring on enhanced mixing. Each measure is asso- 
ciated with a specific physical setting and is amenable to different mathematical 
analysis and/or approximation. Approximations adopting distinct noncommuta- 
tive limiting procedures for control parameters may yield conflicting predictions. 
We then thoroughly investigate a fundamental example of this problem, that of a 
simple steady shear flow stirring a scalar sustained by a monochromatic source-sink 
distribution, to expose an apparent contradiction between the mixing measures 
when the high-Peclet and large length-scale-separation limits are exchanged. Fi- 
nally, we propose a mathematical framework that utilizes information about single 
particle dispersion, a generalization of Batchelor's 1949 theory [1] of diffusion in a 
field of homogeneous isotropic turbulence, to accurately predict scalar concentra- 
tion variance reduction by the (in this case, inhomogeneous and anisotropic) flow. 
This new approach produces a uniformly valid dependence of the effective diffusion 
on the control parameters, reconciling the apparent inconsistencies. 

2. Mixing measures. A conventional modeling approach, and the one we will 

focus on here, is to describe the flow-enhanced mixing by replacing the advection- 
diffusion operator with a effective diffusive operator, i.e., 

kAT-u-VT -> V-(K eff -VT) (3) 

where K eff is an effective diffusivity tensor designed to capture some specific fea- 
ture^) of the mixing process. 

One such feature is transient passive particle dispersion with K eff defined by 

Ktf = \ Um f t n(X t (t) - Xi(0))(X,-(t) - X,(0))], i,j = l,-..,d (4) 

where X(t) = (Xi(t), • • • ,X^(t)) T is the trajectory described by (1) and the sta- 
tistical average E[-] is taken over all Brownian paths W(t). The effective diffusion 
tensor Kij is often defined as lim^oo ^E[(JQ(£) —Xi(0))(Xj(t) —Xj(0))], but these 
are equivalent when the covariance elements converge grow no more than linearly 
with time. 
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The steady sinusoidal shear u = iy/2U s'm(k u y), where i is the unit vector in 
the x— direction, is a case in point. The effective diffusivity tensor defined by the 
long-time dispersive behavior of particles is the Taylor dispersion result [16] 

K e ff= ^d + Pe 2 ) O^j (5) 

where the strength of the advection is gauged by the non-dimensional Peclet number 
Pe = UI u /k with the flow-characteristic length-scale l u = k~ x . 

In the absence of steady sources and sinks, and when the separation between the 
characteristic length-scale of the scalar density distribution, call it Id, and that of 
the flow, l u , is very large, homogenization theory [5, 8, 9] asserts that the scalar 
concentration evolves according to a diffusion equation with effective diffusion tensor 
K eff for a broad collection of deterministic and stochastic flows. The normalized 
tensor ft -1 K eff emerges naturally as a dimensionless measure of the mixing efficacy 
of the stirring, and previous analysis for the ld/l u oo homogenization limit has 
shown that each component in the normalized tensor is bounded above in terms of 
the Peclet number according to 

if eff 

£HT := ^< 1 + p^ (6) 

Hi 

Taylor dispersion, and in particular the steady sinusoidal shear flow, saturates the 
bound for Kff . 

Homogenization theory explicitly assumes the large length-scale separation be- 
tween the flow and the scalar field. In one incarnation a steady large-scale gradient is 
imposed in one spatial direction to formulate the so-called "cell problem" [6]. This 
is implemented theoretically by writing T = —Gx + so the advection-diffusion 
equation (2) becomes the evolution equation for the concentration fluctuation field 
0: 

d t 6 + u-X76 = kA6 + G(i- u). (7) 

The flow and fluctuation fields are then assumed to be periodic and mean-zero on 
a cell of size l u across which the "mean" scalar gradient G is held constant. This 
model serves as the starting point for many studies of turbulent mixing [13]. 

The solution of (7) provides another, potentially distinct, measure of mixing en- 
hancement in terms of the scalar flux-gradient relationship. If the effective diffusion 
coefficient Kff is defined as the ratio of the mean scalar flux in the x— direction to 
the mean scalar gradient in the x— direction, then 

;= Kg_$. ( «r - - vp) _ 1 (jvffi < 1 p<? 

Hi HiG G z 

where (•) denotes the long time and spatial average within a periodic cell. The 
second expression for in terms of follows from the time-cell average of 

times (7) which implies that Hi{\V0\ 2 ) = G(0(i • u)), and the upper bound follows 
from 

(|W| 2 ) _ (0(1 ; u)) (yg- v-^-g)) 

s (IWlWv-'»l¥ =Pe(|V()|a)l , 
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where the Peclet number is defined in terms of the rms velocity U and the charac- 
teristic length-scale of the flow l u as 

U 2 :=(\u\ 2 ), ll:=(\V-^)/(\u\ 2 ). (10) 

Here the inverse gradient operator V -1 has the Fourier symbol — i\k\~ 2 k operating 
on mean-zero functions. The same shear flow u = i\/2U s'm(k u y) that saturates the 
bound in (6) also saturates the bound in (8) because the steady-state solution to (7) 
is 0(y) = \/2GU sin(k u y) / tvk 2 . In homogenization theory, the effective diffusivity 
defined in (4) is often computed by solving (7). 

Another mixing measure was recently introduced motivated by problems with 
spatially inhomogeneous scalar sources and sinks [3, 11, 14, 17]. Transporting 
particles from sources to sinks by advection may help to suppress the inevitable 
inhomogeneities in the scalar field beyond that which molecular diffusion can ac- 
complish alone. Inhomogeneities in the scalar concentration may be measured by 
its space-time averaged variance, and stirring efficacies can be defined in terms of 
the suppression of scalar variance on various spatial length-scales. In particular, in 
the presence of a given source-sink distribution t), Thiffeault et al. introduced 
the notion of scale dependent equivalent diffusivities [3, 17] 



Kp ~ v (iv^i 2 ) ~ K y (iv^p) (n> 

where the unstirred "reference" scalar field 6q satisfies 

d t 6 = kA6 + s (12) 
while the stirred scalar field satisfies 

d t 6 + u- V6 = kA6 + s. (13) 
The associated mixing efficacies 

are measures of the statistical steady-state flow-enhanced concentration variance 
reduction field on small (p = 1), intermediate (p = 0), and large (p = —1) spatial 
scales. Here the source-sink distribution s(x,t) is, without loss of generality, spa- 
tially mean 0. This notion parametrizes the flow-enhanced mixing by the effective 
or "equivalent" molecular diffusion coefficient that achieves the same level of scalar 
concentration variance suppression that the stirring supplies. 

Rigorous analysis shows that given a stationary source-sink distribution s(x), the 
efficacies and equivalent diffusivities are bounded by 

£ p <CPe (15) 

as Pe — )> oo for a wide class of deterministic and stochastic flows with a prefactor 
C depending on l 3 /l u , where l s is a length-scale characterizing the source-sink dis- 
tribution. This linear-in-Pe bound can be saturated for some source-sink and flow 
combinations, as shown by Plasting and Young [11]. Moreover, for time-dependent 
flows with statistical homogeneity and isotropy properties often associated with 
fully developed turbulence, even smaller estimates for £ p , i.e., asymptotic upper 
bounds < Pe a with a < 1, hold for certain classes of source-sink distributions and 
depend on the spatial dimension [3]. These estimates have also been shown to be 
sharp [10]. 
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The discrepancy between the Pe 2 scaling (6) and the Pe 1 scaling (15) of the 
diffusion enhancement factors casts doubt on the universal applicability of the ap- 
proaches and their associated mixing measures, and raises questions about modeling 
mixing by (7) when Pe is large. This discrepancy may to some extent be attributed 
to the necessity of a scale separation between the tracer concentration and flow 
fields in the basic particle dispersion analysis, and its inevitable implication in flux- 
gradient models. Equivalently, maximally enhanced diffusion in these approaches 
requires time to develop, time which may be as large as 0(l 2 u j n\ the effective "mean 
free time" for a typical tracer from the initial distribution to move by molecular 
diffusion onto a streamline in another direction. Said differently, Taylor dispersion 
may require a long time to emerge [18]. 

In the presence of sustained scalar sources and sinks this separation of time scales 
may never effectively be achieved: a relevant time scale in the sourced problem is 
the lesser of Z^/ft, the time for enhanced Taylor-like dispersion to appear, and l s /U , 
the time it takes for a particle to be advected by the flow from a source to a sink. 
If l s /U < Iu/k, the bulk of the scalar variance may be dominated by particles 
that are most recently replenished and depleted rather than by particles that have 
been in the system for a long time, and are thus relatively well mixed. Transient 
features of particle dispersion cannot be neglected and the simple parameterization 
of the advection-diffusion operator with a stationary tensor may not be valid. As a 
result, mixing efficacies for source-sink problems may differ from those deduced from 
transient mixing problems when l s /U < I 2 / ft, i.e., when at high Peclet numbers 
when Pe > l s /l u . The specific example analyzed in detail in the next section 
precisely illustrates this noncommutativity of the large length-scale-separation limit 
and the large Peclet number limit. These two distinct asymptotic limits can produce 
different predictions for the effective diffusion scaling. 

3. Explicit example: Sinusoidal source-sink and shearing. Consider the 
simplest nontrivial case where the steady single length-scale flow 

u(x,t) = iV2U sm(k u y) (16) 

is stirring the steady single-scale source 

s(x,i) = V2Ssin(k s x) (17) 

as depicted in Figure 1. The non-dimensional control parameters are the Peclet 
number Pe and the length-scale separation ratio r, which we define as 

Fe = -^< r = T = r- ( 18 ) 

The reference, unstirred, steady scalar concentration is 0o(x) = sm(k s x) while 
the exact steady (t — >• oo) stirred solution oo (x : y) solving 

y/2 U sm(k u y)d x oo = kAO^ + V2S sm(k s x) (19) 

has the form 

6>oo(x, y) = f(y) sin(k s x) + g(y) cos(k s x) (20) 

where / and g are I 21 — periodic functions that can be computed via numerical or 
asymptotic methods. 

Different dependencies of the different mixing enhancement measures at large 
Peclet numbers are by now well documented [9, 14, 16, 18]. In the limit r — >• oo the 
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Figure 1. The sinusoidal shear u = i\/2U sm{k u y) with source 
distribution s(x,t) = y/2 S sm(k s x) in the shaded background. 




Figure 2. vs k 2 u nt for k u y = 0, f , f . 
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direct approximate solution obtained by homogenization theory, i.e., the solution 
of the steady inhomogeneous diffusion equations effective diffusivity given by the 
particle dispersion or flux-gradient definition, is simply 

V2 S sm(k s x) Oq(x) 
° MX) = ^(1 + Pe 2 ) = T~fP?' (21) 

This suggests that the homogenization theory approximation is 



,HT_ (IV^OI 2 ) 1 | p 2 (22) 

for all p, although careful application of homogenization arguments to the gradient 
of the scalar field reduces the homogenization theory prediction for E\ to O(Pe) as 
Pe^oo [7]. 

To see (22) we derive the explicit formula for tracer particle position covariance 
from the solutions of the stochastic differential equations written as 

X(t) = x + y/2U /J sm(k u Y(s))ds + V^ W ± (t), 
Y(t)=y + V^W 2 (t) 

where X(t = 0) = (xo,yo) T and W\ and W2 are two independent Brownian motions. 
The entries in the 2-by-2 covariance tensor G(yo,t) are computed as 

Cnfoo,*) = E[(X(t)-x ) 2 ] 

= 2a^ + 2/7 2 / E[sin(fe n y'(s))sin(/c n y"(r))]dsdr 

J[0,t}2 

= 2Kt + ^4 - 1 + e _/c ^ £ — cos(2/c n ?/o) 

U 2 

~ 2(«+— J)*, as t^oo, (24) 
nk u 

C 12 (y ,t) = E[(X(t)-x )(Y(t)-y )] 

= 2U^ f E[W 2 (t)sm(k u Y(s))]ds 







C 22 (j/o,t) = E[(F(t)-^ ) 2 ] = (26) 

The homogenization theory approximation (5), i.e, classical Taylor dispersion, fol- 
lows from the long time limits of these quantities. Figure 2 shows the temporal 
evolution of K\\ for different values of yo. Note that it takes time « 100 /k 2 K for 
the full Pe 2 enhancement to emerge. 

In contrast, it has been shown [17] that the efficacy So of any flow on the torus, 
steady or time-dependent, stirring the simple monochromatic source-sink distribu- 
tion is bounded according to 

S < Vl + r 2 Pe 2 . (27) 




Because rPe = U / nk Sl this upper bound does not depend on any length-scales in 
the flow. In fact, for this source-sink distribution this bound is saturated by a 
steady, spatially uniform wind directly blowing source to sink and sink to source. 

For the sinusoidal shear flow u(x, t) = iy/2U sm(k u y), a detailed high-Pe asymp- 
totic analysis for the exact steady-state solution [14], what we will refer to here 
as Internal Layer Theory (ILT), shows 

£0 ~ r 7 / 6 Pe 5/e , S^-rPe, £ +1 ~ r^Pe 1 ' 2 (28) 

as Pe — )> oo at fixed r. Comparing the asymptotic bounds (22) and (28), it is 
clear that with two control parameters, Pe and r, the large- Pe asymptotics does 
not commute with large-r asymptotics. As a result, the dependence of £q on large 
Peclet numbers has two distinguished regimes that cross over near r = Pe. 

Figure 3 illustrates the accuracy of different theoretical estimates for £q by com- 
parison with the exact value for 10 _1 < Pe < 10 6 and r = 10 m , m = —1, • • • ,6. 
The exact behavior of So shows that the two limits, Pe — >• oo and r — >• oo, do not 
generally commute and thus lead to two distinct asymptotic regimes of the Pe de- 
pendence: the 1 + Pe 2 regime for Pe < r and the r 7 / 6 Pe 5 ^ 6 regime for Pe > r. The 
question we now turn to is whether transient particle dispersion information can be 
utilized to correctly predict steady-state concentration variance suppression in the 
presence of steady sources and sinks. 



4. Dispersion-diffusion theory. In order to reconcile the notions of effective 
diffusion in terms of particle dispersion, on the one hand, and source-sink sustained 
scalar variance suppression, on the other hand, we propose what we call Dispersion- 
Diffusion Theory (DDT). Specifically, we retain the dependence of the effective 



MIXING MEASURES AND EFFECTIVE DIFFUSION 



9 



diffusivity tensor on time [1] and initial location [18] by defining 

K(fo,i) := Ijt ^ ^ = ljt n(Xi{t) ~ Xi (°M X i® ~ X ^- ^ 
Then without sources and sinks, it is hypothesized that the probability density of 
a single passive particle may be approximated by the diffusion equation 

d d d 

—p(x,t]x ,t ) = —Kij(xo,t- to)—p(x,t]x ,t ) (30) 

with initial distribution p(x, to] xq, to) = 5(x — xo). For spatially periodic problems 
with period L, this fundamental solution is 

p(x,t',x ,t ) = Yl e^^-^ o) -^- c( ^' £ - £o) ^. (31) 

Now we propose to approximate the solution to (13) with a source-sink distribution 
s(x,t) by the integral 

#ddt(^,£)= / dt / dxop(x,t;xo,t )s(x ,to) . (32) 
Jo J[o,L] d 

That is, we simply apply the principle of linear superposition to the particle den- 
sity introduced (or depleted) at all positions xq at all past times to. Note: it is 
straightforward to show that #ddt does not satisfy an "effective" advection-diffusion 
equation, even if s(x) is steady. 

This approximation models each individual tracer particle's position by a Gauss- 
ian probability distribution with the proper variance (and here, with mean zero, 
although mean displacements could be included as well). Dispersion-Diffusion The- 
ory generalizes Batchelor's 1949 theory [1] for stirring by homogeneous turbulence 
by retaining the initial position dependence in the effective diffusivity; this matters 
for inhomogeneous flows such as the steady sinusoidal shear flow. 

The exponential decay in the kernel (31) suggests that the integral (32) may 
be dominated by the behavior of p for small t — to. This is the mathematical 
implementation of the physical statement that the bulk features of the scalar field 
are determined by the particles most recently injected and depleted by s, although 
this feature is not uniform in the wavenumbers k (i.e., in the relevant length-scales). 

To evaluate the DDT approximation we define 



f(y)+im ~ T dt U e <nfc«(y-yo)-i(fc. l nfc tt )C(yo,t)(fc. l nfc 14 r d2/dB 

(33) 

Then from (31), (32) and (24-26) we have 

<Wr(^, y) = f(y) s'm(k s x) + g(y) cos(k s x) (34) 

with real functions 

f(y) = Yl cos ( nk uy)ln, g(y) = Y sin ( nk uy)ln (35) 

s n— — oo s n— — oo 

and the dimensionless integrals 

I n = h^L r dt cos(nfc„yo)e-^>" fe ") c ^ ^ nk ^ T dy , n e Z. (36) 
27T ,/n l n 
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As an approximation to the exact solution #oo, we use #ddt to compute the 
multiscale mixing measures and to evaluate its ability to recover the parameter 
dependences of the efficacies for large Peclet number and/or large scale separation. 
From (33) through (36) and Parseval's Formula, the approximate multi-scale mixing 
efficacies are 




m \ r Pe 



^-i/^P^- VE jf r2n2)/r (37) 



r 



Pe 



T 11 

m l+r 2 n 2 

where the sums $^ n (-) are taken over all integer values of n. 

In the limit of large scale separation, r = k u /k s — )• oo, a straightforward change 

of variables to z = k u yo and r = u ^ - in (36) using (24-26) suggests 

£„^o=^ (^) 

Thus, the intermediate and large scale mixing efficacies approximated by #ddt 
satisfy 

£ ~ £-1 ~ Pe 2 (39) 

in the limit r — >• oo for large but finite Pe in agreement with the homogenization 
theory prediction 1 + Pe 2 . 

As shown in Figure 4, the dispersion-diffusion approximation to the efficacy 
(DDT) is visually indistinguishable from the exact value. While the homogeniza- 
tion theory prediction (HT, dot-dashed line) applies up to Pe w r and the internal 
layer asymptotic approximation (ILT, dotted line) is accurate for Pe > max{r, 1}, 
the dispersion-diffusion approximation is uniformly accurate. 

Note that the dispersion-diffusion approximation (along with the exact results, 
of course) respect the rigorous efficacy bound (27). The efficacies for a range of r 
are plotted vs. rPe in Figure 5. It is interesting to observe that the simple sine 
flow nearly saturates the absolute upper bound, which holds for all possible stirring 
flows, when r w Pe. We expect the dispersion-diffusion approximation to respect 
the efficacy bound more generally as well. This is because at high Pe we expect the 
major contribution to #ddt(^^) in (32) to come from integration times to within 
l s /U of t. Tracer particle position variance may reasonably be (upper) estimated 

by 

E[(Xi(t) - XMKXjit) - Xj(0))] Z (2«< + U 2 t 2 ) (40) 
as t — >• 0, so for steady sources and sinks each Fourier mode may be estimated 

|4tWI>^ (41) 

as Pe — )> oo. This implies that 

cDDT Ul source I source Ul u //io\ 

o ^ " - j x —— - r ^ e \^ Z ) 

rv in rv 



MIXING MEASURES AND EFFECTIVE DIFFUSION 



11 



10 12 



(MO 



<0 



— - Exact 

ioial — DDT 

10 8 
10 6 
10 4 
10 2 
10 0J — 



id 



id 2 

Pe = 



id 4 



r = 10 6 
r = 10 5 

r = 10 4 

r = 10 3 

r = 100 

r = 10 

r = 1 

r = 0.1 



u 



Figure 4. £ vs Pe, for r = 10 m , m 



10 r = 100_ r 10 3 



10 4 



-1,0,- 
= 10 5 



10 6 



,6. 



10 b 



10 12 




Figure 5. £o and the rigorous upper bound (solid line) vs rPe, 
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where the distinguished length-scale characterizing a general source-sink distribu- 
tion is 

^source ^|y-l s |2^ ' V ' 

as long as it is no n- vanishing. Should I source vanish, i.e., if the source-sink distribu- 
tion contains too many small length-scale components, we would expect different 
(sublinear) high-Pe scaling [3]. 

To evaluate the potential for the various theories to capture the mixing efficacies 
at large and small scales, we plot the different £_i and £+\ and their approximations 
in Figure 6 and 7. For the large-scale mixing efficacy the homogenization 
approximation is accurate up to Pe ~ r while the dispersion-diffusion approximation 
appears to capture the correct scaling for any Pe modulo a constant prefactor error. 
The direct homogenization approximation fails to capture the correct behavior of the 
small-scale efficiency £ + i for any Pe > 1, which is not unexpected since it explicitly 
neglects small scale structure in the scalar field. But a careful homogenization 
analysis focusing on the gradient of the scalar field [7] predicts £+i ~ Pe in the 
1 < Pe < r regime and this does correctly capture the intermediate behavior. 
The DDT approximation, on the other hand, appears to follow the direct (naive) 
homogenization approximation for Pe < r but then adopts the correct scaling when 
Pe ^> r, albeit with a prefactor error. 

To visualize the detailed structures in the scalar field as captured by the various 
approximations, Figure 8 compares the exact and the #ddt and #ht approxi- 
mations for Pe = 10 m , m = 1, • • • , 5 for fixed r = 1000. The panels in each column 
(fixed Pe) are plotted in the same grayscale from —1 (black) to 1 (white) where 
the scalar fields are normalized by the magnitude of the sup-norm of the corre- 
sponding exact solution, ||#oo||oo- As Peclet increases with the intensity of the flow, 
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Figure 7. E\ vs Pe, for r = 10, 100. For large r, I n in (37) is 
very sensitive to numerical errors in evaluating £+\ producing the 
artificial fluctuations seen in the DDT curves. The KKS-HT curve 
refers to the homogenization theory result for the gradient of the 
scalar field by Keating, Kramer, and Smith [7]. 



internal layers develop in the strongly-sheared regions where the speed is relatively 
small; in the weakly-sheared regions, however, the speed is large and the scalar 
field is well mixed and has small variance. It is clear from the figure that although 
the dispersion-diffusion approximation fails to capture the details of the layers, it 
does reveal the bulk behavior of the scalar field and more importantly, recovers 
the correct bulk scalar variance. The homogenization theory approximation greatly 
overestimates the effect of the stirring when Pe> r. 

5. Conclusion. By direct comparison of various mathematical models and mea- 
sures of mixing we have shown that different definitions of effective diffusion predict 
distinct large Pe scalings of flow-enhanced diffusion. The discrepancies result from 
the diverse physical mechanisms motivating the different definitions of the mixing 
measures. In particular, 

• In the transient mixing problem (s = 0) the long-time behavior of the well- 
dispersed scalar density is controlled by the advection and diffusion of tracer 
particles ignoring scalar density structure on length-scales of the stirring. As 
a result, the effect of the stirring may be described by a diffusion equation 
with an effective diffusion tensor enhanced by as much as a factor of Pe 2 . 

• In the presence of sources and sinks, the scalar concentration generally de- 
pends on the length-scales of both the flow and the source-sink distribution. 
We cannot generally use the long-time, large length-scale dispersion results to 
approximate the system with an effective diffusion equation: homogenization 
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Pe=10 Pe=100 Pe=1000 Pe=10000 Pe=100000 




Figure 8. Scalar fields for (x,y) e [0,27r/fc s ] x [0,2ir/k u ] when r = 1000. 



theory is not applicable without the pristine separation of scales between the 
stirring and the source-sink distribution. The enhancement of molecular dif- 
fusion by stirring generally depends in a nontrivial way on both Pe and the 
scale separation r. 

• Dispersion-Diffusion Theory, motivated by the desire to utilize the essential 
information in the particle dispersion process in the presence of sources and 
sinks, reconciles the non-commutative limiting procedures adopted in the lit- 
erature. DDT retains the dynamical and inhomogeneous aspects of the ef- 
fective diffusivity tensor used in homogenization approach and approximates 
the scalar concentration with an integral similar to the solution of a diffusion 
equation with sources and sinks. 

• The Dispersion-Diffusion approximation should generally respect the upper 
bound in (15). Indeed, at high Peclet number the dominant contribution to 
6 ddt comes from the most recent times which, due to the gaussian nature of 
the approximation, leads to variance suppression < Pe 1 . 

We may thus utilize the classical particle dispersion perspective to accurately pre- 
dict enhanced mixing via scalar variance suppression, even for highly anisotropic 
and inhomogeneous flows and without a separation of length-scales. This is crucial 
for problems where the sources and sinks possess the same range of scales as the 
stirring. We do not, however, yet see how to uniformly reconcile the predictions 
of the flux-gradient model in (7). This model is frequently adopted as the defin- 
ing framework for turbulent stirring and mixing, but it incorporates infinite scale 
separation from the start so no such reconciliation may be possible. This raises 
the question of the relevance of the flux-gradient model to applications involving 
statistically steady state mixing in the presence of sources and sinks. 

Dispersion-Diffusion Theory may be applied to more general source-sink distribu- 
tions [3] and/or more complicated flows like homogeneous and isotropic turbulence 
where, following Batchelor [1], the dispersive behavior of passive particles is modeled 
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by 

E[(Xi(t) - Xi^Xjit) - Xj(0))} ~ (2Kt + UH 2 + . . . )Sij , (44) 

at least for displacements within the inert ial range. The term in the covariance 
~ t is due to molecular diffusion while the term ~ t 2 characterizes the short term 
drift. For source-sink distributions with a well-defined spatial scale falling below 
some "outer" scale l u of the turbulence (i.e., r < 1), a calculation very similar to 
that in (40) through (43) with the simple dispersion relation (44) yields 

£ p ~ r Pe, (45) 

for Pe ^> 1 > r, which saturates the scaling of the rigorous upper bounds. This 
suggests that homogeneous isotropic turbulence may be a nearly optimal mixer in 
this sense of steady state scalar variance reduction. An important aspect of the 
future research is to test this conjecture with direct numerical simulations and/or 
experiments for passive scalars that are sustained by steady sources and sinks while 
being stirred by turbulent flows. 
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